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ao.  ABSTRACT  (Continue  an  toootoo  aide  II  necoreary  and  Identity  A y  Alack  nuaAor) 

■  In  this  work  a  promising  new  electron  cyclotron  maser  oscillator  is 
proposed  and  analyzed.  The  configuration  utilizes  an  open  resonator  cavity 
containing  a  gyrating  electron  beam  which  translates  along  an  external 
magnetic  field.  The  magnetic  field  is  directed  perpendicular  to  the  axis 
of  symmetry  of  the  open  resonator.  Because  the  wave-particle  Interaction 
volume  is  extremely  large,  the  total  input  electron  beam  power  can  be  high 


and  the  power  density  low. 
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20.  Abstract  (Continued) 

This  configuration  has  a  number  of  potentially  very  attractive  features. 
Among  them  are:  i)  high  operating  radiation  power  levels  (  >  MWs),  11)  high 
efficiency  operation  (>45%),  ill)  naturally  suited  to  short  wavelength 
operation  (X<2mm),  iv)  operates  efficiently  at  low  electron  beam  voltage 
(10-100  keV),  v)  natural  transverse  mode  selectivity,  vi)  moderate  insen¬ 
sitivity  to  beam  temperature  effects.  The  non-linear  interaction  between 
the  electrons  and  resonator  field  have  been  analyzed  and  an  expression  for 
the  steady  state  efficiency  obtained. 

The  expression  for  efficiency  has  a  rather  simple  analytic  form  when 
the  amplitude  of  the  resonator  fields  are  small,  i.e.,  small  signal  regime. 

In  this  case,  the  efficiency  is  an  essentially  odd,  non-oscillating  function 
of  the  frequency  mismatch  (difference  between  resonator  frequency  and 
relativistic  cyclotron  frequency). 

For  a  uniform  external  magnetic  field,  total  efficiencies  in  excess  of 
30%  can  be  realized.  We  have  considered  enhancing  the  efficiencies  by 
spatially  contouring  the  magnetic  field.  Appropriately  contouring  the 
magnetic  field  across  the  resonator  by  ~  5%  Increases  the  efficiency  from 
~  35%  to  ~45%. 
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theory  of  the  quasi -optical  electron  cyclotron  maser 

I.  INTRODUCTION 

At  the  present  time,  two  classical  radiation  mechanisms,  i.e.,  the 
free  electron  laser (FEL)  and  electron  cyclotron  maser (ECM),are 
under  extensive  study  because  of  the  great  potential  they  show  as  new 
classes  of  coherent  radiation  sources.  Experimental  results  on  the 
pel (20  22)  gQ^(22  25)  ^avfi  veTy  encouraging. 

In  free  electron  lasers  the  active  medium  is  a  beam  of  relativistic 
electrons.  Such  sources  have  the  potential  for  generating  coherent  radia¬ 
tion  ranging  from  the  millimeter  to  the  optical  regime  and  beyond.  They 
are  frequency  tunable  and  in  principle  extremely  efficient  generators  of 
intense  radiation. 

The  electron  cyclotron  maser  in  its  present  form  has  reached  a  far 
more  mature  stage  of  development  than  the  FEL.  In  the  millimeter  regime, 
electron  cyclotron  masers  have  generated  power  levels  substantially  higher 
and  more  efficiently  than  the  more  conventional  radiation  sources.  Experi¬ 
mental  efficiencies  are  impressive,  e.g. ,  22%  efficiency  at  x  *  2mm  with  a 
CW  output  power  of  22  kW. 

In  this  paper  we  propose  and  analyze  a  new  electron  cyclotron  maser 
oscillator  configuration,  which  utilizes  an  open  resonator  cavity.  Our 
quasi-optical  cyclotron  maser  has  a  unique  potential  for  becoming  a  new 
type  of  coherent  radiation  source.  In  principle,  the  device  is  capable 
of  generating  coherent  radiation  in  the  millimeter  to  stfcmillimeter  regime. 

Manuscript  submitted  October  6, 1980. 
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at  power  levels  in  excess  of  magawatts,  with  efficiencies  exceeding  SOI. 

The  basic  structure  consists  of  an  open  resonator  containing  a  beam  of 
electrons  gyrating  about,  as  well  as  streaming  parallel  to,  an  applied 
magnetic  field.  The  magnetic  field  is  directed  transverse  to  the  axis  of 
the  open  resonator,  which  consists  of  two  or  more  appropriately  curved 
mirrors.  Moderately  low  electron  beam  energies  can  be  used,  i.e.,  10-100 
keV,  even  though  the  wave-particle  interaction  mechanism  is  due  to  relativ¬ 
istic  effects.  This  configuration  has  a  number  of  distinct  advantages  over 
the  more  conventional  radiation  sources.  Some  of  these  advantages  are: 
i)  extremely  high  operating  power  levels,  ii)  high  operating  frequency, 
iii)  high  efficiency  and,  iv)  natural  transverse  mode  selection.  Since  we 
utilize  an  open  resonator  and,  thus,  have  a  large  interaction  voltme,  the 
input  electron  beam  power  can  be  extremely  high  while  the  power  density  can 
be  kept  moderately  low.  The  usual  limitations  on  beam  power  imposed  by 
space  charge  effects  can  therefore  be  overcome.  Since  the  wave -particle 
interaction  is  fairly  efficient  (~50%) ,  high  radiation  power  levels  can  be 
achieved.  The  operating  frequency  is  limited  solely  by  the  external 
magnetic  field  and  is  independent  of  the  dimensions  of  any  physical 
structure.  Favorable  coupling  between  the  electrons  and  radiation  field 
occurs  near  harmonics  of  the  relativistic  cyclotron  frequency. 

A  quasi-optical  resonator  has  many  modes  which  in  principle  can  experi¬ 
ence  gain,  producing  a  multi-mode  output  signal.  The  fundamental  transverse 
resonator  mode  can  be  preferentially  excited  in  the  open  resonator.  If 
the  mirrors  in  the  open  resonator  are  made  large  enough  to  intercept  a 
large  fraction  of  the  flux  in  the  fundamental  mode,  we  may  expect  this 
mode  to  have  a  large  Q.  The  higher  order  modes  can  be  expected  to  have  a 
substantially  smaller  Q  since  they  suffer  from  diffraction  losses.  This 
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is  an  advantage  over  a  closed  cavity,  where  all  modes  would  have  comparable 
values  of  Q.  Longitudinal  modes  in  resonance  within  the  resonator  can 
undergo  gain  and  produce  a  multi-mode  output  signal,  unless  they  are 
suppressed.  Longitudinal  mode  selection  can  be  achieved  by  employing  a 
mode  selection;  one  such  selector  is  the  Smith-Fox  interferometer.^^ 

This  new  maser  configuration,  like  the  conventional  electron  cyclotron 
maser  (gyrotron),  has  a  wide  range  of  practical  applications.  These  appli¬ 
cations  range  from  electron  cyclotron  heating  of  fusion  plasmas  to  advanced 
radar  and  communications  systems.  Because  of  the  high  field  level,  short 
wavelength  and  extended  interaction  volume,  the  quasi-optical  maser  may  be 
a  natural  electro-magnetic  pump  source  for  a  free  electron  laser.  In  this 
application  a  second  highly  relativistic  electron  beam  propagating  along  the 
axis  of  the  open  resonator  would  interact  with  the  resonator  field  and 
induce  high  frequency  radiation.  The  frequency  of  the  scattered  radiation 
would  be  ~  4-yy2u) ,  where  is  the  longitudinal  gamma  factor  of  the  second 
relativistic  electron  beam  and  u>  is  the  frequency  of  the  quasi-optical 
maser. 
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II.  NON-LINEAR  ANALYSIS  OF  QUASI-OPTICAL  MASER 


The  quasi -optical  cyclotron  maser  configuration  is  depicted  in  Figure  1. 
The  lowest  order  mode  is  the  well  known  TEMQQ  Gaussian  radiation  bean^  see 
Figure  2.  For  the  electric  field  vector  primarily  polarized  in  the  x 
direction,  the  field  components  of  this  mode  are 

Ex(x,y,z,t)  =  E  (x,y,z)  sin  (k^y  +  ot(x,y,z))  cos  wt,  (la) 

Bz(s,y,z,t)  =  E  (x,y,z)  cos  (kyy  +  a(x,y,z))  sin  a.t,  (lb) 

where  ky=a./c,  E(x,y,z)=E0(rQ/rs(y))exp(-(x2+z2)/r2s(y)) ,  ot(x,y,z)= 

R_1(y)  (x2+z2)  (i»)/c)/2  -  tan^Cy/y^, 

w  is  the  radiation  frequency,  EQ  is  the  field  amplitude  at  the  origin,  rQ 

is  the  minimum  spot  size  at  the  plane  y=0,  r$(y)  =  rQ  (l+y'/yR  )  is 

2 

the  spot  size  at  the  plane  y,  yR  *  rQ u>/2c  is  the  Rayleigh  length  and 

2  2 

R(y)  =  y(l+yR  /y  )  is  the  radius  of  curvature  of  the  spherical  wavefront 

at  y.  The  y  components  of  the  field  are  E  =  -  (c/w)9  E  /  ax  and 

y  * 

By=  -  (c/a.)  3BZ/  9z.  Note  that  |Ey|,|By|  «  |EX|,  |BZ|.  One  can  show  that 
the  E  and  B  fields  satisfy  the  appropriate  boundary  conditions  on  the 
mirror.  This  is  equivalent  to  saying  that  the  mirror  radius  of  curvature 
equals  the  wavefront  radius  of  curvature. 

The  intermode  frequency  spacing  is  6u.  *  ttc/L  where  L  is  the 
separation  between  the  mirrors  and  the  diffraction  angle  of  the  radiation 
is  0d  =  */*rQ  where  A  -  2  ttc/u.  is  the  radiation  wavelength. 

Other  orientations  of  the  external  magnetic  field  and,  hence,  the 
gyrating  electric  beam  are  possible.  Careful  analysis,  however,  shows 
that  these  other  configurations  may  not  be  as  straightforward  to  inclement. 


For  example,  another  possible  configuration  is  where  the  magnetic  field 
and  the  streaming  electrons  are  directed  along  the  y  axis.  It  can  be 
shown  that  to  maintain  wave-particle  coherence  the  thermal  spread  in  the 
y  component  of  electron  velocity  must  satisfy  5vy/vy  «  X/2ttL.  For  long 
resonators  this  condition  places  a  rather  stringent  requirement  on  the 
electron  beam  quality. 

In  Figure  1,  the  gyrating  beam  electrons  rotate  in  the  x-y  plane 
and  stream  along  the  external  magnetic  field  BQ  which  is  directed  along 
the  7.  axis.  For  convenience  we  locate  the  sheet  electron  beam  of  width 
Lb  on  the  y-z  plane, as  shown  in  Figure  1.  Furthermore,  for  reasons  of 
analytic  simp! icity.we  take  the  guiding  centers  of  the  electrons,  upon 
entering  the  resonator  fields,  to  lie  on  the  y-z  plane.  Also  the  electrons, 
upon  entering  the  resonator,  are  assumed  to  have  the  same  transverse  and 
parallel  velocities. 

Strong  coupling  between  the  electrons  and  resonator  field  will  occur 
at  frequencies  near  multiples  of  the  relativistic  electron  cyclotron  fre¬ 
quency.  Let  us  consider  the  fundamental  cyclotron  interaction, 
where  fiQ=  |e|B0/m0c,  y  =  (l+p*p/!UpC2  )**  and  p  is  the  electron  momentum 
vector.  The  electron  Lanror  radius  is  in  general  much  less  than  the 
radiation  wavelength,  i.e.,  r^=  YVj./^0-  t  «  X,  where  vx  =  cB  ,  is 

the  transverse  electron  velocity.  The  minimum  radiation  spot  size  is 
much  greater  than  the  Lanror  radius,  r0>^>>rL*  Therefore,  by  choosing 
the  width  of  the  electron  beam  to  be  somewhat  less  than  the  Rayleigh 
length,  <  i^,  the  resonator  fields  in  (1)  felt  by  the  electrons  can 
be  accurately  approximated  by  the  plane  wave  fields. 
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Ex(y,z,t)  =  E(z)  sin  (|y)  cos  (wt)  , 
By(y,z,t)  =  E(z)  cos  (|y)  sin  (u>t)  > 

where  E  (z)  =  EQ  exp  (-z2/rj  )  and  Ey  =  By=  o. 


We  now  express  the  particle  momenta  g  and  transverse  position  (x,y) 
as  functions  of  Lagrangian  independent  variables.  In  general,  a  convenient 
set  of  Lagrangian  independent  variables  for  this  problem  are  the  z  position 
of  the  particles,  the  initial  momentum  space  angle  0Q,  initial  transverse 
coordinates  of  the  guiding  center  xgQ  and  ygQ  and  the  particle  entrance 
time  into  the  resonator  field  tQ.  Since  the  resonator  fields  fall  off 
like  exp(-z2/rQ2),  the  entrance  position  of  the  sheet  beam  z^n  can  be 
taken  to  be  a  few  spot  sizes  away  from  the  y  axis,  lzjnl>>r0-  In  our 
present  analysis,  xgQ=  0,  |ygQ|  <  Lb  and  tQ  is  the  time  the  particle 
crosses  the  z  =  zin  plane.  The  functional  dependence  of  the  particle 
momenta  vector  and  transverse  coordinates  is  p  =  i?(z »ygo*0o»to) , 
x  «  x(z,yg0,00,t0)  and  y  =  y(z,yg0,©0,t0) . 

The  orbit  equations  for  the  electrons  are 

Pz  dpy  -  -  | e  j  (ym0Ex+p  (Bq+Bz)/c),  (3a) 

dz 

Pz  &  •  1*1  Px  <VBz>/c  •  (3I>) 

£■  0.  (3C) 

where  the  fields  Ex  and  Bz  are  given  by  (2),  with  t  replaced  by  the 

Lagrangian  time  variable  t  (z,  ygQ,  0Q,  tQ)  *  tQ  +  A  /vz  and  vz= 

2in 

pz(z’ygo,0o,to)/Y(z,ygo’Go’to)  mo  is  the  1on9ltudinal  particle  velocity. 


Note  that  p,  is  a  constant  of  the  motion  denoted  by  p,.,  hence,  v, f  0 
unless  pZQ  =  0.  The  Lagrangian  time  variable,  t,  is  the  time  it  takes 
a  particle  to  arrive  at  z  if  it  crossed  the  z=z.jn  plane  at  time  tQ  with 
a  momentum  space  angle  0Q  and  guiding  center  position  (*gO=0,  ygO=0).  Since 
we  are  considering  only  the  fundamental  cyclotron  interaction,  an  appro¬ 
priate  representation  of  the  solutions  of  the  orbit  equations  in  (3)  are 

Px(z)  =  Pxg(z)  +  Pj2)  cos  (wt+0),  (4a) 

P  (z)  =  pyg(z)  +  pA(z)  sin  (uvr+0),  (4b) 

where  PXg»Pyg  are  the  components  of  momenta  associated  with  the  guiding 

centers,  p^is  the  transverse  particle  momentum  and  wt+0  is  the  particle 

momentum  space  angle.  In  (4)  the  dependent  variables  Pxg>Pyg>  Px»  and 

0  are  assumed  to  be  slowly  varying  functions  of  z  as  well  as  functions  of 

ygQ,  0q  and  t^.  By  "slowly  varying"  we  mean  that  the  quantity  has  no  high 

(cyclotron)  frequency  Fourier  components.  The  variables  p  ,p  ,  p  and 

yp 

O  are  not  functionally  independent  of  each  other.  In  fact,  by  requiring 
that  they  be  slowly  varying  functions  of  z,  we  will  derive  four  separate 
but  coupled  equations  which  uniquely  determine  them.  The  field  amplitude 
E(z)  defined  in  Eq  (2),  denotes  the  profile  of  the  radiation  beam  and  is 
a  slowly  varying  function  of  z  since  the  electrons  undergo  many  cyclotron 
orbits  while  traversing  the  resonator  fields.  The  initial  values  of 
dependent  variables  upon  entering  the  open  resonator  fields  are 

‘W2-zin>  ■  Pyg(pzin>  '  °Jz'zfn>  *  Kv  ,Z(z"z1n>*VP»>/Y°mo  and 
0(z=z.jn)  =  0Q.  Furthermore,  it  will  be  shown  that  the  guiding  center 

drift  momentum  of  the  particles  is  much  less  than  the  transverse  momentum, 

i.e.,  I Pxg I  » I Pyg I  <<  PA •  Noting  the  form  for  py  in  (4b),  we  see  that  the 
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(5) 


dependent  variable  y  can  be  approximated  by 

Pj. 


cos  (arr+0)  , 


0  0 


where  yg  is  the  slowly  varying  y  position  of  the  guiding  center,  hence, 

Pyg  =  PZQdyg/dz  and  =  rL  1S  the  electron  Laraior  radius. 

We  now  substitute  (4a, b)  together  with  (5)  into  the  orbit  equations 
(3a, b)  and  carry  out  the  appropriate  operations. Keeping  in  mind  that  E(z) 
is  a  slowly  varying  function  of  z  we  equate  rapidly  varying  terms  and 
slowly  varying  terms  and  discard  all  frequencies  of  order  2f2Q.  These  manip¬ 
ulations  yield  four  interrelated  self-consistent  non-linear  equations  describing 

the  spatial  evolution  of  p  ,  0,  p  and  p  ,  When  r, «A ,  these  equations  are 

xg  yg  l 


dp^=  - 1 e | Eym  ["  sin  (|  y  )  cos  0+cos  (^y  )(^f  cos  0  +  sin  0)] 
di  ^7°  [  c  9  c  g  mQc  mQc  J 


J[6a) 


d|  =  ( mo/pzo  +  ^i^o  sin  (f  yg)  sin  0  -  cos  (£  yg) 


sin  0  -  COS  G) 
'me  me 

o  o 


*1  b 

* 


p 

xg 


ryg 


te'E  p-*  cos  (£  y  )  sin  G 

2W 


lEp.»  <  1  -  ™  >  cos  (|y  )  cos  0  , 


(6b) 

(6c) 

(6d) 


Tfi  m„c 
o  o 


where  dyq/dz  =  Pyg/Pz0  and  Y=(1+(PZ0  +  pj1  J/n^c2)*5.  Equations  (6)  com¬ 
pletely  describe  the  non-linear  steady  state  particle  dynamics  for  the 
fundamental  cyclotron  interaction.  Since  |E/B0(«l,it  is  noted  from  (6c) 


and  (6d)  that  I Pxg I » I Pyg t<<:  Pj.*  justifying  the  corresponding  earlier 
assumption.  The  trajectory  of  each  electron  is  described  by  the  set  of 
equations  in  (6).  The  initial  conditions,  however,  are  different  for  each 
electron,  as  required  by  the  Input  distribution.  For  an  entering  cold  un¬ 
bounded  electron  beam,  the  initial  conditions  are  such  that  at  z=z^n, 

pxg=0  *  pyg=0,  p-l=  pz=pzo*  ^o  where  Qo  ran9es  from  0  t0  2"  and 

yg=ygo  where  ygo  ranges  from  “V2  t0  Lb/2‘  while  pjo  and  pzo  are  the 
same  for  all  electrons. 

The  equations  in  (6)  can  be  considerably  simplified  by  noting  from 

(6c,d)  that  Pxg/m0c  »  0(YVaE/2BQ)  «1  and  Pyg/moc%0  (P  W^cw).  Hence, 

the  second  and  third  terms  in  the  brackets  on  the  right  hand  side  of 

(6a, b)  can  be  neglected  and  y  in  (6a ,b)  can  be  replaced  by  y  .  The 

g  ■'  go 

resulting  equations  are 


These  equations  are  very  similar  to  those  analyzed  in  our  initial  studies 
of  the  non-linear  behavior  of  the  cyclotron  maser  instability. 

In  a  temporal  steady  state  oscillator  the  efficiency  of  converting 
beam  power  to  radiation  power  Is  given  by 

n  '  (Pb,iI.-pb,out>/pb,in  («) 

where  ^  and  Pf,  out  is  the  total  electron  beam  power  flowing  into  and 
out  of  the  open  resonator  in  the  z  direction.  For  the  cold  beam 
distribution  presented  earlier,  the  efficiency,  as  defined  in  (8),  can 
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be  shown  to  give 


n  = 


-  e 


VV^zo1" 


Using  (2a)  and  (4a)  we  find  that 
E  p 


Z  =  z0(jt  27t/u)  Lb/2  2tt 

^  f dz  /  2I7S  fi°ExPx-  {9) 

-Lb/2  6  * 


ExPx  = 


sin  (^  ygo)  cos  0  , 


(10) 


where  we  have  used  the  approximations  y«ygQ  and  Pxg  =  0.  Substituting  (10) 
into  (9)  and  noting  that  Expx  is  independent  of  tQ  we  obtain  for  the 

efficiency  the  expression 


_  -|e|  a) _ 

4tty„(y„-1  )v, 

0  0  zo  0 


/ 


z=z 

.  out 

4/2 

dz 

T 

zin 

-V2 


/ 


2tt 

dO 

I 

2tt 


0  E  (z)  £(z.Vygo)si"(c  VC0S  e(z,eo,ygo). 


(11) 


The  expression  in  (11)  gives  the  full  non-linear  steady  state  operating 
efficiency.  Before  solving  (11)  in  the  fully  non-linear  regime  it  is 
illuminating  to  first  solve  the  orbit  equations  in  the  linear  regime 
and  thereby  obtain  the  analytic  form  for  the  efficiency.  The  orbit 
equation  in  (7)  can  be  linearized  by  setting  pA  =  pi0^?!1^  and  0=0^+0^] 
where  pj0^  and  0^  are  zero  order  quantities  in  the  field  amplitude,  E, 
and  Pj.^  and  o^are  first  order  quantities.  Solving  (7)  we  find  that 


Pi  =  Px0  . 
(0) 

0  (z)  =  0 


(z-zin). 


(12a) 

(12b) 
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(12c) 


z 


(o) 

COS  0 


(i)  lei 


±o  sin(-y  )  / dz  /dz  E(z)cos  0 


zin  zin 


//)+  1-  9^  (12d) 

p  90 

jo  o 


where  Aw  =  w  -  ^0/y0  is  the  frequency  mismatch. 

Now  substituting  (12)  into  (11)  and  carrying  out  the  ygQ  integra¬ 
tion  we  arrive  at  the  following  expression  for  the  linear  efficiency 


=  -(|e|/m0 


„2\2  /*ZOUt  /*Z  ( 

oj^  ydz  y<Jz  I  E(z)E(z')cos(|  (z- 


|  (z-z)  ^a-) 
c  6zo 


w/c  /dz  E(z)E(z#)sin(|  (z-z)  ^j/“)  [  ,  (13) 

V  C  S“  ’ 

in 

where  E(z)  denotes  the  profile  of  the  radiation  beam.  In  the  present 

2  2 

case  the  radiation  beam  has  a  r-aussian  profile  given  by  E(z)=E0  exp(-z/r0). 
Substituting  this  form  for  E (z )  into  (13)  and  taking  the  beam  entrance 
and  exit  planes  to  be  respectively  z^n=-  *  and  zQut=  <*>,  the  integrations 
over  z,z  and  z  can  be  carried  out  analytically.  The  overall  expression 
for  the  linear,  small  signal,  efficiency  takes  the  rather  simple  form 


7T  Y0  /Of 

"'^V1  VV 


-Aw/w)2 

«o  * 


[%  $  *  ->] 


where  =  (r0w/c)/820.  Unlike  the  case  of  a  beam  propagating  along  the 
z  axis,  the  structure  of  n  is  non-oscillating  in,say,rQ.  The  reason  for 
this  is  that  the  radiation  field  is  a  smooth  function  of  z  and  has  no 


abrupt  change  in  behavior  at  the  entry  and  exit  points  of  the  electrons. 
From  (14)  we  see  that  the  efficiency  is  positive  when  &a)/aj>2/(Bj.0C0)2 , 
hence,  the  output  frequency  is  always  slightly  higher  than  the  relativistic 
cyclotron  frequency.  For  typical  choices  of  parameters,  the  term 
B20£^(Aw/(jj)/2  is  much  greater  than  unity.  In  this  case  the  linear 
efficiency  maximizes  when  the  frequency  shift  is  Au)  =  w/C0. 

Appendix  A  contains  a  derivation  of  the  small-signal  efficiency 
at  the  fundamental  as  well  as  at  all  cyclotron  harmonics  using  the 
linearized  Vlasov  equation.  In  this  appendix  the  linear  efficiency 
expression  is  derived  for  the  cases  where  the  electric  field  of  the 
radiation  beam  is  polarized  in  the  x  direction,  i.e.,  polarization 
considered  in  the  body  of  this  paper,  as  well  as  in  the  z  direction. 

The  former  polarization  of  the  electric  field  is  shown  to  result  in 
substantially  higher  linear  efficiencies  compared  to  the  latter  polari¬ 
zation.  In  conventional  electron  cyclotron  masers  the  particle  inter¬ 
acts  with  the  TE  mode  of  the  structure.  This  corresponds,  in  our 
present  configuration,  to  the  electric  field  polarized  in  the  x 
direction.  Polarization  of  the  electric  field  in  the  z  direction 
would  correspond  to  a  TM  mode  interaction  in  the  conventional  configuration. 

In  order  for  the  system  to  operate  in  the  assumed  steady  state, 
there  must  be  losses  which  just  compensate  for  the  power  loss  of  the 
electron  beam.  These  losses  are  composed  of  both  the  output  radiated 
power,  and  the  real  losses  due  to  diffraction  or  dissipation  in  the 
mirrors.  All  of  these  losses  are  characterized  by  the  Q  of  the  cavity 
so  that  the  total  power  out  of  the  resonator  is  pout=w ^stored^’  w^ere 
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^stored  is  the  energy  stored  in  the  cavity,  which  is  proportional  to  E*. 
Since  linear  theory  shows  that  power  lost  by  the  beam  is  also  proportional 
to  Eq,  the  operating  field  amplitude  cannot  be  determined  by  linear  theory 
alone.  Thus  the  oscillator,  unlike  the  amplifier,  is  an  inherently  non¬ 
linear  device. 


As  will  be  shown  in  the  next  section,  where  the  nonlinear  electron  dy¬ 
namics  are  calculated, the  power  lost  by  the  beam  ultimately  levels  off  as  E* 
increases.  Thus  the  actual  operating  point  can  be  calculated  from  the  inter¬ 


section  of  the  graphs  of  power  loss  by  the  beam,  and  power  lost  by  the 
cavity,  as  functions  of  E*,  shown  schematically  in  Figure  3.  It  is 
apparent  that  steady  state  can  be  achieved  only  if  the  losses  are  small 
enough  that  the  two  curves  intersect.  Also,  for  the  optical  cavity 
configuration,  we  find  that  the  power  lost  by  the  oeam  never  becomes 
negative  as  E*  increases.  Thus  steady  state  operation  is  not  possible 
without  power  loss  by  the  cavity. 

The  threshold  condition  for  starting  the  oscillations  in  the 
resonator  is 

n  Pb,in  -  u^tored^  * 


(15) 


where^tored  *  ^<E*/8tt)  irr*L  is  the  stored  field  energy.  To  obtain  the 
threshold  electron  beam  power,  necessary  to  start  the  resonator,  we  use 
the  small-signal  efficiency  in  Eq.  (14).  Substituting  the  maximum  small- 
signal  efficiency,  i.e.,  when  C0Wu=l,  Into  (15)  we  find  that  the 
product  of  the  beam  power  and  resonator  Q  needed  to  start  the  oscilla¬ 
tions  is  PDfinQ  >  4.6x10*  Y0(Y0-l)6j0/62i0[watts]  .  (16) 
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Throughout  our  analysis  on  the  quasi -optical  maser  Me  have  made 
the  tacit  assumption  that  the  energy  lost  by  the  electrons  goes  Into 
supporting  the  assumed  Gaussian  radiation  beam.  This  assumption  is 
common  to  all  conventional  oscillator  problems  and  has  been  proven 
valid  experimentally.  Though  we  have  not  rigorously  proven  this  point, 
concerning  our  present  configuration,  we  have  assumed  it  applies  here  too. 
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III.  NUMERICAL  RESULTS  AND  ILLUSTRATIONS 

a)  Efficiency  of  Maser  in  Uniform  External  Magnetic  Field 

As  an  illustration  of  the  potential  operation  of  the  quasi- 
optical  maser  we  choose  a  60  keV  (yQ  *  1.118)  electron  beam  having  initial 
velocity  components  =  0.4  and  &ZQ=  0.2.  The  location  and  curvature 
of  the  mirrors  are  chosen  so  that  the  radiation  spot  size  is  5.9  wave¬ 
lengths,  i.e.  rQ=  5.9x.  The  linear  (small  signal)  efficiency  given  by 
(14)  is  shown  in  Figure  (4)  for  various  values  of  the  normalized  field 
amplitude,  E0/B0-  The  linear  efficiency  has  a  single  positive  maximum 
when  CqAw/ui  *  1  which  corresponds  to  the  frequency  U)max=(l_  1/5o>'Vv 
The  intermode  frequency  separation  of  a  resonator  of  length  L  is 
6(d  =  ttc/L.  Therefore,  a  spectrum  of  natural  modes  can  exist  within 
the  resonator  separated  in  frequency  by  6id.  When  the  oscillator  is 
first  started  up,  the  mode  frequency  closest  to  uj^ax  will  be  excited 
first  and  grow  to  a  level  where  it  suppresses  the  slower  growing 
natural  modes  further  away  from  u^ax.  It  turns  out  that  the  maximum 
non-linear  efficiency  occurs  at  a  frequency  slightly  higher  than  u^ax, 
i.e.,  for  C0Aw/ui  somewhat  larger  than  unity. 

In  the  absence  of  longitudinal  mode  selection  or  equivalent 
scheme, it  would  not  be  possible  to  take  advantage  of  the  higher  non¬ 
linear  efficiencies  occurring  at  frequencies  greater  than  u^ax  To 
suppress  these  unwanted  modes  near  <i>max  and  operate  the  maser  at  the 

frequency  of  maximum  non-linear  efficiency,  we  will  assume  that  a  longi- 

26 

tudlnal  mode  selector  Is  employed.  Another  approach  seems  possible: 
we  can  start  the  oscillator  at  the  natural  operating  frequency  u^ax. 

When  the  mode  saturates,  the  external  magnetic  field  can  be  decreased 
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slightly  so  that  S0Aw/ui  >  l  has  the  appropriate  value  to  maximize  the 
non-linear  efficiency.  We  will  assume  that,  by  using  either  of  these 
approaches,  the  operating  frequency  can  be  freely  chosen. 

The  non-linear  efficiency  defined  in  Eq.  (11),  is  evaluated  by 
following  numerically  the  particle  trajectories  according  to  (6).  For 
sufficiently  small  field  amplitudes,  the  simulations  reproduced  very 
accurately  the  linear  efficiency  in  (13).  The  beam  entry  and  exit  points 
were  taken  at  +2rQ  and  the  electron  orbits  were  integrated  using  a  4- 
point  Runge-Kutta  integrator.  Figure  (5)  demonstrates  the  higher  non¬ 
linear  efficiencies  achievable  at  larger  values  of  ZqAw/u.  The  efficiency 
in  this  figure  is  obtained  by  solving  (11)  as  a  function  of  EQ/B0  for 
various  values  £oWu>.  A  maximum  efficiency  of  33%  is  obtained  for 
C0Au)/a)  =  7  at  E0/Bq  =  2.25  x  10"2.  Figures  (6)  and  (7)  show  the  spatial 
variation  of  the  efficiency  within  the  resonator  for  various  values  of 
the  normalized  frequency  shift  and  normalized  field  amplitude 

Eq/Bq.  The  spatial  oscillations  in  efficiency  within  the  resonator 
are  due  to  the  oscillations  of  the  trapped  particle  distribution.  The 
Gaussian  profile  of  the  radiation  beam  is  also  shown  in  Figure  (5)  for 
reference  purposes. 

The  examples  presented  so  far  should  not  lead  to  the  impression 

that  good  performance  is  necessarily  associated  with  frequency  mismatches 

significantly  larger  than  the  values  corresponding  to  maximum  linear  gain. 

For  exaiqple,  for  a  beam  with  p,|Q  ■  0.1  and  *  0.2,  corresponding  to 

Y0  *  1.026  (13.3  keV),  the  nonlinear  efficiency  is  plotted  in  Fig.  8  against 

the  radiation  field  amplitude  for  various  values  of  the  ratio  r  A.  In  all 

o 

cases  the  frequency  mismatch  was  taken  to  correspond  to  the  naxinun  linear 
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gain  value,  obtained  from  Eq.  (14).  As  can  be  seen  in  Fig.  8,  a  peak  effi¬ 
ciency  of  20%  is  obtained  for  rQA  ■  2.5,  while  the  values  of  n  ■  17.51  are 
obtained  for  rQA  =1.6  and  3.2.  For  these  cases,  the  frequency  mismatch  is 
given  by  Auu/u  =  0.75  x  10'^,  1.30  x  10"^  and  0.58  x  10’^,  respectively,  and 
corresponds  to  optimal  linear  gain  value. 

b)  Efficiency  Enhancement  By  Contouring  External  Magnetic  Field 

It  is  possible  to  enhance  the  non-linear  operating  efficiency  of 
the  quasi-optical  maser  by  either  pre- bunching  the  electron  beam  in  momentum 
phase  angle  uur  -0  or  by  appropriately  contouring  the  external  magnetic  field. 
Pre-bunching  the  electron  beam,  by  utilizing  a  two  open  resonator  Klystron 
type  configuration, is  in  principle  straight forward  and  results  in  extremely 
high  efficiencies.  However,  depending  on  the  length  of  the  ballistic  phase 
bunching  region  (distance  between  the  two  resonators) ,  electron  beam  thermal 
effects  may  present  a  problem  at  high  frequencies.  Contouring  the  magnetic 
field  appears  to  be  the  simplest  method  for  enhancing  efficiency.  By  slightly 
contouring  the  magnetic  field,  as  a  function  of  z,  a  more  advantageous  momentum 
phase  distribution  of  the  electrons  can  be  realized  with  a  single  resonator. 

A  significant  improvement  in  efficiency  over  the  already  highly  efficient 
uniform  magnetic  field  case  can  be  realized  in  this  way.  Figure  9  shows  the 
spatial  evolution  of  efficiency  with  and  without  magnetic  field  contouring. 

The  magnetic  field  in  this  case  was  decreased  linearly  by  5%  between  the  points 
z  -  -2rQ  and  z  =  2tq.  For  this  variation  a  final  total  efficiency  of  45%  was 
achieved. 

c)  Design  Examples 

We  conclude  with  two  specific  detailed  design  exanples ,  which  will 
demonstrate  the  potential  of  the  configuration  we  have  analyzed.  In  the 


17 


following  examples  the  steady  state  performance  of  the  maser,  operating  at 

150  (Hz  (X  =  2  mm),  is  analyzed.  The  electron  beam  is  taken  to  be  generated 

2 

by  a  diode  with  a  current  density  of  10  A/cm  .  Within  the  oscillator,  the 

27 

electron  beam  has  a  rectangular  cross  section,  with  area  A  =  2rQy^,  extending 

from  x  =  -r  /2  to  x  =  +rQ/2  and  from  y  =  -y^  to  y  =  +y^,  where  ro  is  the  spot 

size  and  y^  is  the  Rayleigh  length. 

The  first  exanple  deals  with  the  60  keV  beam  with  =0.4  and 

0  =  0.2,  for  which  the  efficiency  is  shown  in  Fig.  5.  For  tq/\  =  5.9,  the 

2 

beam  cross  section  is  A  =  52  cm  ,  giving  a  beam  current  of  520  A  and  an  input 

beam  power  of  31  MW.  Due  to  the  slight  variations  of  the  radiation  field 

amplitude  across  the  beam,  the  conversion  efficiency  has  to  be  appropriately 

,  -2 

averaged.  For  a  radiation  field  of  magnitude  E  /B  =  2.4  x  10  at  the  center 

of  the  beam,  this  weighted  average  efficiency  is  n  =  28%,  hence  the  radiated 

output  power  is  Prad  *  8.7  MW.  For  such  performance,  a  normalized  frequency 

mismatch  =  7  is  required,  hence  &d/ju  =  3.77%  and  uu/Qq  =  0.93,  For 

X  =  2rnn(u)=9.4x  1011  sec'1),  a  magnetic  field  of  B  =58  kG  is  required. 

o 

In  addition,  by  appropriately  tapering  the  external  magnetic  field  the  output 
power  can  be  increased  to  Pra{j  =  12  NW.  Finally,  if  either  of  the  above 
illustrations  (with  or  without  the  external  field  taper)  is  to  be  achieved 
by  initially  adjusting  the  external  magnetic  field  to  the  value  corresponding 
to  maximum  linear  gain  of  the  operating  frequency,  then  for  a  44  cm  long 
resonator  (i.e. ,  equal  to  2  y^) ,  the  value  of  Q  obtained  from  (16)  must  exceed 
the  180,  while  the  operating  value  of  Q  is  40,000. 

In  the  above  example  the  frequency  mismatch  was  substantially  higher 
than  the  value  required  for  maximum  linear  gain.  The  highly  impressive 
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performance  associated  with  the  adaptation  of  such  a  condition  warrants  its 
implementation.  The  realization  of  such  a  large  frequency  mismatch  would 
require  a  mode  selector  or  an  external  source  to  set  up  the  radiation  field. 
However,  if  such  complications  are  to  be  avoided,  excellent  performance  can 
be  achieved  as  will  be  shown.  In  the  second  example  we  consider  a  13.3  keV 
beam  with  g  =  0.1  and  giQ  *0.2,  interacting  with  a  radiation  beam  of 
X  =  0.2  cm  and  rQA  =  2.5,  at  the  frequency  mismatch  associated  with  maximum 

linear  gain,  i.e.  Aau/u  =  0.75%.  In  this  case,  the  cross-sectional  area  of  the 

2  2 
beam  is  A  *  4  cm  .  Assuming  that  the  diode  current  density  of  10  A/cm  can 

2 

be  conpressed  to  50 A/cm  ,  the  input  beam  power  is  2. 7 MW.  For 
Eq/Bo  *  2.3  x  10~\  the  average  efficiency  is  18%,  hence  the  radiation  power 
is  Pra(j  =  0.50**.  The  required  value  of  the  external  magnetic  field  is 
Bq  *  54.5  kG.  For  this  case,  the  linear  threshold  condition  (16)  is  not 
restrictive,  since  it  sinply  requires  that  the  nonlinear  efficiency  be 
smaller  than  the  linear  value,  which  is  the  case.  Assuning  an  oscillator 
length  equal  to  16  yR  *  63  cm,  the  operating  value  Q  is  1560,  larger  than 
the  threshold  value  of  150,  required  by  (16). 
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APPENDIX  A 


Here  we  work  out  the  steady  state  linear  theory  for  the  energy 
loss  of  the  electron  beam  as  it  traverses  the  wave  fields  in  the  open 
resonator.  In  steady  state,  the  energy  equation  for  the  beam  is 

A  wf  =  E  •  J  ,  (/ 

where  is  the  beam  energy  flux  and  J  is  the  beam  current.  The  con¬ 
figuration  is  as  shown  in  Figure  1.  The  change  in  beam  power  as  it 
crosses  the  resonator  then  is  Jd3r  E*J,  so  the  main  problem  is  to  calcu¬ 
late  an  expression  for  the  perturbed  beam  current  density  J. 

This  can  most  easily  be  done  by  analyzing  the  linearized  Vlasov 


equation 


•lu  +  v-v  - 


'  leIP 


where  a  time  dependence  e  is  assumed.  The  quantity  f^0-*  is  the 
unperturbed  distribution  function,  that  is,  the  distribution  function 
at  z=-“  ,  and  f^  is  the  perturbation  to  it  induced  by  the  fields  in 
the  resonator.  Instead  of  using  independent  variables  r,  £,  it  is  more 
convenient  to  use  as  independent  variables  the  quantities  (xg,yg,z,pA,<j>,pz) 


where 


Px  =  Px  cos  (^o  z  +  <f>) 

Py  *  P*  sin  (^o%  z  +  $)  , 

x  -  x  +  PAsin(2o%  z  +  0) 
x  xg  pz _ 


y  s  yg  .  PxCOS(^0  z  +  0) 


■>  v> 
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where  is  the  nonrelativistic  cyclotron  frequency.  These  have  the 
advantage  of  being  constant  in  the  absence  of  the  radiation  fields  in 
the  resonator.  In  these  new  variables,  the  Vlasov  equation  reduces  to 


(- 4'"'- 4  •  # 


3f 

3p 


(o) 


(„)  ( ° ) 

and  assuming  that  f'o;  is  independent  of  d> ,  ^  reduces  to 


3f 

3p 


(o). 


(2^c»s(?f» 2  **) ifc +  n^iir  £-)*  *  *  *) 


3 


ft :  m  3x_ 

oo  g 


— <L_\  +  e  -i- 
3xg  )  z  3Pz 


r  (  0  ) 


There  are  two  possible  polarizations  for  the  radiation.  The 
first  is 

6E  =  E(z)  cos  ky  e"iu)t  ex+c.c. 


6B 


_  ike  E(z)  sin  ky  e'1wt  ez+c.c.  > 

U) 


which  we  call  TE  since  it  has  6E  A  BQ.  The  second  polarization 

<5E  =  E(z)  cos  ky  e"1a)t  ez+c.c. 

<5B  ="i^-  E(z)  sin  ky  e"ia)t  e  +c.c.  » 

-  (0  X 

is  TM  since  6B  ±  b^. 

For  TE  polarization,  the  linearized  Vlasov  equation  reduces  to 


The  perturbed  distribution  function  then  is 

/  \  /  rz  -iym  wz' 

f(l'  =  exP  (i  “p - j  J  dz'  exp  — p -  G(z'), 

-OO 

The  quantity  we  are  interested  in  is  #SE*‘6J  d3r  which  is 

/■'Q  m  z  '  . 

P.d;\dpzd0dxgdygdz  (p/va^cosC-—-  +  <J>) 


(A10) 


In  doing  this  integral,  exp(iky),  which  appears  in  the  expression  for  6E 

must  be  written  as  exp{  ik  [y  -  (p  /ft  m  )  cos  (ft  m  z/p  +  which  is 

g  J-  OO  OO  z 

equal  to 

?  °n  fe) e1"’  ’(kV  n1[?  ^  *  *])• 

Then  /  E* • J  d3r  is  given  by 

J  oo 

Ji’r  E*-J  -  \ ydxgdygdz  p^dM+dP,  E*(z>^JJn'(vio)exp 

fz'  e«P  [i  ^-a]  E(z)^Jn  (v^)exp^  )] 


n,  /ftjri  z 


Pz  cosl  P2  + 


*  *)(2^  -Hr) 


+  (k-+-k)  +  c.c.  , 


(All) 


where  we  have  assumed  that  ~ —  =  0,  that  is, the  distribution  of  guiding 

% 

centers  is  uniform  in  the  y  direction  (actually  the  dimension  of  the  beam 
in  the  y  direction  is  very  long  compared  to  a  wavelength),  and  also 
have  exploited  the  fact  that  no  quantity  except  f^  depends  on  x  ,  so 
if,  0  term  integrates  to  zero.  Since  f^  is  independent  of  4>,  the 

3*g 


only  term  which  does  not  integrate  to  zero  over  <j>  is  the  zero  Fourier 
harmonic.  This  collapses  the  double  summation  over  n  and  n'into  a 


single  summation, so  that 


E**0=2tt  |e|^xgdygdz  pxdp.,dpz  E*(z)^l 


Jdz  '  exp ji 


4 1 

P-* 


■ 

o^m,. 


n= 

2 


p,  (z-z')fP 


-.  fit-  II 


3p 


fn(v^] 

o) 

+  c.c. 


(A12) 


where  the  (k-*-k)  has  now  been  explicitly  included. 

The  next  problem  is  to  do  the  zz' integral .  Assuming  E(z)  has  the 
form  Eq  exp  (-zVr2),  this  integral  is  of  the  form 

J&/.  Jdz'  expj-  j^r0j  +  +  i  a(z-z ')j| where  a=  (ym^  -  n  mQao)/pz. 

-  oo  -oo 

By  completing  the  square,  the  z'integral  can  be  done  in  terms  of  error 
functions  so  it  reduces  to  a  single  integral 


-  00 


Since  the  error  function  is  an  odd  function  of  u,  the  integral  of  the 
term  in  the  souare  brackets  containing  the  error  function  is 


/ 


i  sin(2arQu)e 


erf  u 


which  is  purely  imaginary  and  sums  to  zero  upon  adding 
conjugate.  Thus  only  the  unity  In  the  square  brackets 
the  integral  and  the  total  result  is 


tt  2  -ot2r,2 

7  ro  e  T“° 


the  complex 
contributes  to 
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Hence  the  total  result  is 


This  is  the  final  result  for  the  power  loss  of  the  beam  as  it  traverses 
the  fields  in  the  open  resonator.  To  progress  further,  assume  the  dis¬ 
tribution  function 

f(,)  -  sr  «<el  -  Pi"  >  {Ip2-  i>zo>  {<V  •  (AW 

where  ois  the  surface  charge  density  of  the  sheet  beam.  The  incident 
power  is 

=  Lo(y-I)  pz  c  2  .  (A15) 

f  Y 

Assuming  n=l,  w  *  JL/ y,  and  0.  (x)=  x/2  (i.e.  ,  kp  «n«) 

v)  1  X  0  O 

the  integrals  can  easily  be  done  and  Eq.  (14)  can  be  recovered  for  the 
efficiency.  For  n=2  one  can  also  show  that  the  efficiency  has  the  same 
basic  form  as  Eq.  (14),  but  it  is  multiplied  by  an  overall  factor  of 

!  k2Pl 
*  ^  ’ 

which  is  much  less  than  unity.  Thus  an  interaction  at  the  second 
harmonic  exists,  but  it  is  weaker. 


For  the  TM  mode  an  analogous  calculation  gives  the  result 
/d3r-  E*-J  (  =  -|e|  v/dxgdygdp]dpz^^rj  exp J- Jo  |Y  moU)*nmono 

J  ~  TM  n^»  I  |_  Pz  J 

E°2  J"2  (§v)(?  (‘  '  v)st"*9if  *  ' 


3  8 


(A16) 
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At  each  harmonic,  the  power  loss  of  the  beam  is  smaller  by  a  factor  of 


order 


(%) 


from  what  it  was  for  the  same  harmonic  with  TE  polarization. 


Thus  the  coupling  of  the  TM  mode  with  the  radiation  in  the  optical  cavity 


is  much  weaker  than  for  the  TE  mode. 
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OUTPUT  RADIATION 


Fig.  3  —  Determination  of  the  point  of  nonlinear  steady  state  operation. 
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